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I. INTRODUCTION 

Nonequilibrium dynamics in closed quantum systems, and in particular quantum quenches, have attracted much 
experimental^"— and theoreticalSri^ attention in recent years. There is a growing consensus that integrable models 
exhibit important differences in behaviour as compared to non-integrable onea^. In particular, by now there is ample 
evidence that the stationary state after a quantum quench in an integrable theory is described by a generalized Gibbs 
ensemble (GGE)^ with density matrix 



Pgge 



cxp 




Here H^^^ is the Hamiltonian and iJ^'-' are local^i integrals of motion 



(1.1) 



(1.2) 



The Lagrange multiplies A; are fixed by the requirement that the expectation values of the integrals of motion are 
time-indcpcndcnt 



lim 



= lim 



Tr [pGGEg^'^] 



(1.3) 



Here L is the size of the system under consideration. In practice it is often useful to work with a truncated GGE^, 
where only the y "most local" conservation laws are retained 



(y) 

PtGGE 



1 



Z, 



iv) 

tGGE 



■cxp -£a[^)hW 



1=1 



(1.4) 



Here the a|^^ are fixed by 



lim 



(^o|H(')|^o) 
L 



Tr 



lim 



L 



1 = 1, 



,y- 



(1.5) 



The full GGE is then recovered in the limit y — > oo, after the thermodynamic limit has been taken first. Assuming 
that a given integrable system indeed approaches a stationary state late after a quantum quench, which is described 
by a generalized Gibbs ensemble, important questions are how to construct the GGE in practice, and how to then 
determine expectation values of local operators. It is these questions wc aim to address for the particular case of the 
spin-1/2 Heisenberg XXZ chain. A priori there are four steps: 



1. Determine the local conservation laws. 



2. Calculate their expectation values in the initial state after the quench. 

3. Construct the GGE density matrix in such a way that equations (|1.3p are fulfilled. 

4. Determine the expectation values of local operators in this ensemble. 
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In the following wc address these in turn. As local conservation laws we consider the minimal set obtained from 
the logarithmic derivative of the transfer matrix at the "shift-point"—. It has been recently found that the XXZ 
Hamiltonian in general has local conservation laws that are not obtained in this wa.y^. In principle they could be 
accommodated in our construction as well. However, in order to keep things simple, we restrict our analysis to the 
antiferromagnetically ordered regime of the Heisenberg chain, where, as far as we know, the minimal set of local 
conservation laws is complete. With regards to step 2, we focus on a class of simple quenches, for which the initial 
states are unentangled. We show how to treat these cases analytically. Our method generalizes to weakly entangled 
initial states of matrix product form, but the analysis becomes much more complicated. The GGE density matrix 
is constructed by the quantum transfer matrix method^. The most difficult issue here is what values the Lagrange 
multipliers \j take. We argue that it is possible to completely specify the quantum transfer matrix, without having 
to explicitly calculate the Xj. Finally, GGE expectation values of local operators can be calculated by borrowing the 
results of the Wuppertal group for finite temperature correlators^^i^. 



II. LOCAL INTEGRALS OF MOTION 



Wc consider the XXZ Hamiltonian 

L 



= i E ^I'^hi + + A(<T|a|+i - 1) , (2.1) 

1=1 

where L is even, cr" are Pauli matrices = cr") and we parametrize the anisotropy as 

A = cos7. (2.2) 

It is well known that (j2.ip is solvable by the algebraic Bethe Ansatz method^. Local conservation laws if ('^^ can then 
be obtained from the logarithmic derivative of the transfer matrix t(A) 

H^'^=^(^§^\o^r[^-,X) . (2.3) 

V 7 OA/ A=0 

By definition the conservation laws commute with one another 

=0. (2.4) 
The transfer matrix is constructed by Algebraic Bethe Ansatz and takes the form 
t(z + A) = Tr[Li(A)Li_i(A)...Li(A)], 

^.(A) ^ i±f^+ smh(^) sinh(.7) , 

' 2 sinh(i7+^) 2 sinh(z7+^)^ ' ' ^ ' 

where t^^v^^ are Pauli matrices acting on the auxiliary space, and the trace is taken over the latter. In the following 
we denote indices in the auxiliary and quatum spaces by Roman (a, b) and Greek (a, /3) letters respectively. 



III. EXPECTATION VALUES OF LOCAL INTEGRALS OF MOTION IN THE INITIAL STATE 

Given an intial state |^'o), we aim to determine the expectation values 

{^olH^^^l^o) . (3.1) 
It is convenient to work with the generating function 

n..(A, = i + + A)|*o, = (-2_)' ^ , ,«) 



k= 



where the right hand side follows from (|2.3p . In order to evaluate f2ip,-,(A) we use that, when viewed as a power series 
in A for large L, we have formally 

-(^ + A) - r(z) exp ( -^ ^ f -2-) ' ^i/w) . (3.3) 



k=l 
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This suggests that for large L we have 



T-l(^ + A) = [r(^ + A)]^ 



(3.4) 



in the sense that the power-series expansions in A coincide. These observations lead to the following (approximate) 
expression for the inverse 



\i + X) ~ Tr[A/L(A)A/L_i(A)...Afi(A)], 



M,(A) 



1+T'V| sinh(2^) l-T^'^l sinh(-i7*) 



sinh(-i7* + ^) 2 sinh(-i7* + ^) 



(3.5) 
(3.6) 



The generating function p.2p can then be expressed as 



n*„(A) 



ld_ 

L dx 

L dx 



^{^o\t{i + x)t~\i + \)\^o) 
Sp(^'o|V^t(x,A)...\/i(x,A)|*o) , 



(3.7) 



where Vn{x, A) arc 4x4 matrices with entries (Ki(a;, A))°^ that arc operators acting on the two-dimensional quantum 
space on site n 



{vux,x))z = {L„ix)riAux)r 



(3.8) 



In this notation Vl{x, A) . . . Vi{x, A) is a regular product of 4 x 4 matrices and Sp denotes the usual trace for 4 x 4 
matrices. Let us now assume that |4'o) is a product state 



Then f7ii,„ can be written as 



f^*o(A) 



l*o) 



)_d_ 

L dx 



-1 



Sp 



Y[u,ix,X) 



where 



U,{x,X) = {^i'^\V,{x,X)\¥,'^). 
We now discuss how to implement the above programme for some explicit examples. 



(3.9) 



(3.10) 



(3.11) 



A. Quench from a;,t) 



Our first example is the product state 



^/2 ■ 



(3.12) 



This corresponds to all spins pointing in the x-dircction. This initial state corresponds to a quantum quench in the 
XXZ-chain with an applied transverse magnetic field 



. L L 
^e^) = I E '^^^^'h-i + '^i'^i+i + A(a|a|+i - 1) - ^ E '^i • 



(3.13) 



Preparing the system in the ground state of iJ(oo) gives the initial state (j3.12|) . and the quench is to the integrable 
zero-field Hamiltonian H{0). Using translational invariance we have 



1 d 



x — X 



Sp 



{u{x,x)y 



(3.14) 
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Denoting the largest eigenvalue of U{x, A) by fimaxix, A), this gives 



19 d 

^x,tW ^ (/^max(a;, A)) =— /iniax(a;, A). 

L ox x=\ ox x=\ 



(3.15) 



In the last step we have used that /Xmax(A, A) = 1. The matrix U{x, A) is readily calculated and its largest eigenvalue, 
which for small x, A is very close to 1, is calculated using Mathematica. This results in 



-ij sin(7) 



2 + 2 cosh(7A) + 4 cos(7) 
Matching the power-series expansion around A = with p.2p gives 

(a;,t liy^'^^la;,!) I - A d''~^ 1 + A 



lim 



TV 4 fe*^-! 

The results for H'^^\ iJ^^) and iJ^^) agree with Ref. [Hj. 



^=Ocosh^(Vl - A'^x/2) + A 



(3.16) 



(3.17) 



B. Quench from the Neel state 



We next consider the Neel state 



I Neel) 



t)2,-l ® I i)2j 



(3.18) 



This would be the initial state for a quantum quench from the ground state of an initial XXZ Hamiltonian with 
A = +00 to a final Hamiltonian with finite A. Using translational invariance by two sites, our expression (|3.10p for 
the generating function becomes 



i^Nccl(A) 



]_d_ 

L dx 



Sp 



(C/i(x,A)C/2(a;,A))^/2 



(3.19) 



where 



C/i(x,A) = i(t|Fi(a;,A)|t)i , 
C/2(x,A) = 2a|V^2(x,A)|i)2. 



(3.20) 



These two 4x4 matrices and the largest eigenvalue of Ui{x, X)U2{x, A) (for small x, A) are readily calculated using 
Mathematica 



f^N6cl(A) 
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sin(27) 



2 2cosh(7A) - 1 - cos(27) 
Matching the power series expansion around A = with p.3p then gives 

(Neel|i?('='|Neel) A d''~^ 1 - A^ 



lim 



L 



2 dx 



fe-i 



c=ocosh(Vl - A^x) - A2 



(3.21) 



(3.22) 



C. More general initial states 

In principle our method can accommodate more complicated initial states of matrix-product form 

|*o>=Trtej^^iA,], 



(3.23) 



where Aj is an to x to matrix with entries that are quantum states on site j, and Tr is the trace over the TO-dimensional 
matrix space. Considering such states is however beyond the scope of the present paper. 
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IV. GENERALIZED GIBBS ENSEMBLE AND QUANTUM TRANSFER MATRIX 

Combining (|1.4p with (|2.3p , we can express the density matrix of the truncated generahzed Gibbs ensembles as 

(v) 1 -ly^ , \,(si^-$-y iogT(t+y)\ 

PiGGE = -(^e ^^1.^0. (4.1) 

Agge 

This density matrix can be analyzed by the quantum transfer matrix approach^. Following the analysis of Kliimper 
and Sakai for computing the thermal conductivity at finite temperature in Ref. (49j . we introduce inhomogeneities in 
the transfer matrix and define the ensemble 

_ lim^^_^oo(T~^(^)T(^ + 2ui.n/7)) ■ ■ ■ (t~^(z)t(z + 2un:n/i)) „x 

P{tll;jv,...,tljv;jv} ^ r7 ■ l^-^j 

As transfer matrices with different spectral parameters commute [r(Ai), t(A2)] = 0, we have 

WmAT . pE" 1 l0gT(i+2Ml;N/7)-l0gT(i) 

iimAf->oo K 

P{ui-M ,---,nN;N} ~ y ■ V^-'^) 

In order to achieve 

P{iii,jv,.. = PtGGEi 

(4.4) 

we need to choose the inhomogeneities Mj;7v such that 

j}^]ll^^'^^S'r{i + '2ui,N/l) -logT{i) ^ -i^Xjl^^^—"^ logT(i + y)|^^P . (4.5) 

t=i j=i ^ y 

A sufficient condition for Eq. (|4.5p to hold is that the spectral parameters Uj-N satisfy 

iT..E[/(^ + 2.,W7)-/»] =-^E^.(^d^)'/(^ + ?/)| ^ (4.6) 

for any function /(y) analytic at y = i. We do not have to require further properties if we can find a solution such 
that 

lim Uj.,N = . (4.7) 
By series expanding / about zero, Eq. (14.61) can be rewritten as 

N 



J = l 

A solution of (|4.8p satisfying (|4.7p is given by 



lim > ui,^ = { ^ ^ ' - ^ (4.8) 



Umy+r,N^u\%^Y.'"n^^:^ J-l,---,?/ m = 0, . . . , [iV/j/J - 1 , (4.9) 

Tl=l 

provided that 

E -i^--i^ = -<^)^!A,. (4.10) 

ni,...,„^=i 

Having represented the density matrix of our truncated GGE in terms of an inhomogeneous transfer matrix, we 
wish to generalize the quantum transfer matrix approach for finite temperature correlation functions (soj in order to 
determine expectation values of local operators in the (truncated) GGE. One can show that 

^GGEPGGE = lim Trj jAt^™ {Q) ■ ■ ■ T'f'^Q^ , (4.11) 



6 



where the monodromy matrix T^'^^'^ is given by 



N-l 



(A{^x/2) B{^x/2)\ _ -pr 



N-l-i 



2u. 



i-N 



1 



(4.12) 



and the matrices L and M are defined in Eqs (|2.5[)p.6[) . In the antiferromagnetic regime A > 1 it is customary to 
define A = cosh 77, corresponding to 77 = 17, and change variables as fohows: 



X = —A . 
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yQTA/ ^2^/7) acts as an upper triangular matrix on the vector 



|0> 



(g) • • • (g) 



and wc have 



with 



C(A)|0)=0, A(A) |0) = a(A) |0) 



D{\) |0) = d{\) |0) 



a(A) 



d{\) 



r ^ 

[n 



sinh(A - v^^) 



.^J{ sinh(A - u^y^ - ??) 
sinh A ^ ^ 



.sinh(A + rf), 

The Bethe ansatz equations for the eigenvalues of the quantum transfer matrix take the form 



a{'Wj) 



M 

n 



sinh(wj — Wk + I]) 



d{wj) sinh(wj - Wk - v) 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



We then introduce the auxiliary function 



N 



= TT sinh(A - Vk + rj) 

a(A) sinh(A - i;^ - ?;) 



(4.18) 



where wi, . . . , uat is the solution of the Bethc ansatz equations corresponding to the largest eigenvalue of the quantum 
transfer matrix t^'^^'{x) = Tr[r<3^^^(.T)]. For A > 1 the function l+a(A) has simple zeros at A = Vj (on the imaginary 

axis) and y poles of order N/y at A = u^^^j^ inside the rectangle Q defined by |ReA| < 77/2 and |ImA| < 7r/2. Using 



these analytic properties we obtain the following integral equation for a(A) 



loga(A) 



lim — 'S^ log 

JV-i-oo y 



sinh (A + 77 - sinh(A) 



sinh(A — u^^Ij) sinh(A + 77) 



sinh(277) log(l + a{uj)) 



Q 2711 sinli(A — uj + rj) sinli(A — a; — 77) 



(4.19) 



where C is a rectangular contour with edges parallel to the real axis at ±i7r/2 and to the imaginary axis at ±7 where 
< 7 < ^. We may use Eq. (|4.6p to replace the inhomogeneities by the Lagrange multipliers Xj specifying the 
truncated GGE 



N 



lim — Vlog 



3inh(A + 77 — u.^jv) sinh(A) 



sinh(A — Uj-^^) sinh(A + 77) 



F 



. sinh 77 d 



where 



F[x] 



2 dX 



sinh 77 



sinh A sinh(A + rj) J ^ 



9 + 



(4.20) 



(4.21) 



J=0 
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Finally we have 

logfl(A) = F 



.sinhr; d 



2 dA 



sinhr; 



dw 



sinh(277) log(l + a(u;)) 



sinh A sinli(A + rj) J Jq 27ri sinli(A — w + 77) sinli(A — lo — rf) 



(4.22) 



We stress that this equation is valid only if the expectation value of S*^ is zero. For non-vanising (5'o|5^|5'o) one 
needs to add a constant to the driving term, which corresponds to the Lagrange multiplier of the total spin in the z 
direction (see e.g. Ref. (soj ) 



logfl(A) = h + F 



' sinh 7/ d 
i — 



dA 



sinh 7] 



sinh A sinh (A + ry) 



duj 



sinh(277) log(l + a{uj)) 



Q 2'Ki sinh(A — cj + r/) sinh(A — cj — 77) 



(4.23) 



For the sake of simplicity we will restrict our analysis to the h ^ case in the following and return to the ft, 7^ in 
section I Villi 

The inverse function fl(A) = l/fl(A) fulfils the following integral equation 



loga(A) = F 



.sinh 77 d 



dA 



sinh T] 



sinh(2?/) log(l + a{uj)) 



h(A) sinh(A — rj) ' Jq 2'Ki sinh(A — cj + 77) sinh(A ~ uj — rj) 



(4.24) 



Similarly to the thermal case, the Lagrange multipliers enter into the calculation of correlation functions only through 
the auxiliary functions fl(A), fl(A). 



A. Thermodynamic Properties 



Thermodynamic properties are completely determined by the largest eigenvalue Ao(0) of the quantum transfer 
matrix. The logarithm of Ao(0) is given by^ 



logAo(O)^ hm hm 



Au) sinh77log(l +fl(w)) 
Q 27ri sinh u sinh(aj + 77) 



dw sinh77log(l +a(w)) 
27ri sinh sinh(aj — 77) 



The expectation values of the local conservation laws can then be expressed in the form 



91ogAo(0) 



dX 



iv) 



do; 



sinh 77 dx-a(uj) 



c 2'Ki sinhw sinh(aj + 7;)(1 + fl(a;)) 



dw 



sinh77(5A. loga(a;)) 



Using (|4.24p . we obtain 

„ 1 ^.sinh77y-i 
9 („) loga(A) - 



sinh^ 77 



9 (y) logfl(A) 



2 sinh(A) sinh(A — 77) 

:pres 
sinh 7/ \ i 



dw 



Q 2'Ki sinh w sinh(tt' — 7;)(l+fl(a;)) 
sinh(277)(9^(y) loga(a;)) 



(4.25) 



(4.26) 



Jq 2'Ki sinh(A — w + 77) sinh(A — lo — ri){\ + a(a;)) 

(4.27) 

The solution to this integral equation is conveniently expressed in terms of an auxiliary function G 



di-'GiX,fi;0) 



M=0 



(4.28) 



where 



G{X,^l■,a) = -coth(A-7i) + e"''coth(A-M-??)+ </> |^^^^^^^4^i^(A-cj;a) 



c 2Ki 1 + a(a;) 



K{X;a) = e"" coth(A - 77) - coth(A + 77) 
Substituting this back into (|4.26l) . we obtain 



.f-ismhi]y^j_^ 



sinh 77 G{uj, /i; 0) 



Q 2Ki sinhw sinh(tj — 7;)(1 + a(w)) 



(4.29) 



(4.30) 



For the case of the full generalized Gibbs ensemble, i.e. the limit y ^ 00 oi the truncated GGE considered above, 
we may lift the relationship (|4.30p to the level of the generating function ft^^ (/i) defined in p.2p 



do; 



sinh 77 G{lu, /i; 0) 



(J 27r7 sinh sinh(aj — 77)(1 + a(aj)) V 77 



, 



i = i,2,.. 



(4.31) 
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This can be rewritten in a more compact form as 

do; sinh?7 0(0;, 0) 



c 27ri sinhcj sinli(aj — 77)(1 + a(w)) 77 ° \ V 



(4.32) 



In passing we note that the magnetization can be computed analogously by taking the derivative with respect to 
the magnetic field of the integral equation (j4.23p . which gives 



duj G(w,0;0) 
c 27ri 1 + a{u}) 



(4.33) 



To summarize the result of this section: in the thermodynamic limit the GGE can be formulated in terms of a 
quantum transfer matrix, whose largest eigenvalue is given in terms of the functions a and a, which in turn are defined 
through the nonlinear integral equations (j4.22p and (|4.24p respectively. The expectation values of the local integrals of 
motion can then be expressed in terms of a and the auxiliary function G'(A, /x; 0) defined in (j4.29p . In order to proceed 
for the case of a truncated GGE, it is necessary to solve our system of integral equations subject to the constraints 
()4.30p . Carrying out such a computation entails calculating the Lagrange multipliers A^^'. On the other hand, we 
are ultimately interested in the full GGE itself. In this case it is possible to avoid having to determine the Lagrange 
multipliers A^-^\ as we will show next. 



V. ELIMINATING THE LAGRANGE MULTIPLIERS 

In order to proceed, it is convenient to switch to the "6b-formulation" of the integral equations by defining 



b{x) = afix'^ + ^j , h{x) ^ a(ix 



2, 

9^i^) ~ ±G(ix^±^,ifj.;Oj , x^=x±ie. 



The functions (j5.1[) can be shown to obey the integral equations 

"sinr; d 



\ogh{x) = F 
logb(x) = F 



2 dx 
sin 77 d 



2 dx 



gtix) = ~d{x- n) + 



5p (x) = -d{x - ^l) + 

r'r/2 



d{x) + [k * log(l + h)]{x) - [k^ * log(l + b)](a;) 
d{x) + [k * log(l + h)]{x) - [k+ * log(l + h)]{x) 

(x) 



k * 




at 1 


{x)- 






9^, 1 


1 


+ b-iJ 


* — 

1 


+ b-^J 


k * 




.9; 1 


[x)- 






1 


1 




* — 

1 


+ b-iJ 



where [A * h]{x) = ^h{x - y)h{y) and 

'^^''^ " ^ cosh(?7n) ' ^^""^ " e2„|„| + 1 



k±{x) = k[x^ =b irf) , 



n— — 00 



n— — 00 



(5.1) 



(5.2) 



(5.3) 



As shown in Refs [43, d[x) and k{x) can be expressed in terms of special functions. In terms of the new variables, 
equation (j4.32p is rewritten as 



2 da; 



—d{x) [ + __i ) = 4fc^) + — (-2^/77) . 

TT Vl+b \x) l+b V 



(5.4) 



The Lagrange multipliers Aj are just parameters of the integral equations, which are fixed by Eq. (|5.4p . Quite 
surprisingly, we can remove the explicit dependence on Xj by taking the difference between the first two equations of 
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system (|5.2p . i.e. considering the system of integral equations 



logb(x) -logb(a;) ^ [(fc+ + A:) * log(l + b)](a;) - [(fc- + fc) * log(l + b)](a;) 



Ai. 



4fc(A^) + -f^*o(-W'?) 



-d{x - m) + 
—d{x) 



k * 



1 + 



1+b 



{x)- 



k+ * 

9n {X) 



9, 



A' 



1+b 



l+b" 



(^) , 



.l + b-^(x) l + r'(a;) 

In general (|5.5p have to be solved numerically by iteration. In addition, from the the same equations it follows 

e-""[logb]„ + e""[logb]„ 



where 



i^[m sinh rj] 



[f]r 



2 r\T 



(5.5) 



(5.6) 



(5.7) 



By setting n = in (j5.6p . one obtains the inverse temperature (the Lagrange multiplier of the Hamiltonian) . The 
other Lagrange multipliers are more difficult to compute, as it is not always easy to identify the Taylor coefficients 
of F[x] from the values of the Fourier coefficients on the right hand side (the case considered in Section IVI AI is 
an example). This suggests that solving the system (|5.5p is less demanding than working out (|5.2p by recursively 
computing the Lagrange multipliers. 



VI. PERTURBATIVE EXPANSION 



We now consider examples, in which (j5.5|) can be solved through an expansion that can be carried out analytically. 
Importantly, this shows that it is indeed possible to completely specify the quantum transfer matrix describing the 
GGE without having to calculate the Lagrange multipliers \j. 

If we consider the expectation value on the ground state of the model, it is known that 



(6.1) 



where /c(/i) is given in (|5.3p . Eqn (|6.ip can be inferred directly from (|5.5p . (|5.6p by considering the case b{x) = 
e(l + 0(e^)) and b{x) = e(l + 0(e^)) in the limit e ^ 0. This is equivalent to the zero temperature limit in the Gibbs 
ensemble, which in the A > 1 case clearly tends to the thermodynamic ground state of the model (see also Ref. [50|). 

This observation opens the door for carrying out a perturbative expansion of the system (|5.5p in the limit of a 
"small" quench. We define functions p{x) and ({x) by 



1 



l+b-^{x) 
1 



- pC(a;)/2 



p{x) 



- p-C(^)/2 



1 + b (x) 

In terms of these new variables, (15.51) are rewritten as 



p{x) 



(6.2) 



C(a-) = -[(fc + k+) * log(l - e«/V)](a;) + [{k + fc_) * log(l - e~^/^p)]ix) + log(^ 



p{x) 



-C(-)/2p(x) 



9^{x) = -d{x -p)+\k* (pe«/2g+)l (a,) 



9f, (x) = -d{x - /i) + k* {pe '^/^^^ ) (x) - 



'^d{x)p{x) 



fc_*(pe ^/^g^)J(x) 
k+ * (pe'^/2.g+)l (x) 



--fJa-o (-2^/77) - fc(/i) 
r? 



(6.3) 
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The definition of a "small" quench is one for which \p{x)\ ^ 1. In this case we may solve (|6.3p by iteration. At lowest 
order we have 

p{x)^p^'\x), C{x.)^0, g^{x)^-d{x~f,). (6.4) 
At this order, the last equation of (j5.5|) reads 

[d * (p(i) d)] i^i) ^ 2k{pi) + (-2/./r?) . (6.5) 

This can be solved by Fourier techniques. The n'th Fourier coefficient is 

- ' ■ 2* 17*„(-2;./,), (6.6) 



where 



cosh(77n) l + e^''l"l 77 7_x tt 
and going back to real space we obtain 

(1),, ^»^o(-^+^) + ^l^^o(-^-») + X,,(x) 

j^,(x)= (6.8) 
^ coshT?- cos(2x) ^ ' 

In the next order of iteration one replaces p{x) with p'^^\x) in the equations defining Q{x) and g^(a;), and solving the 
resulting system gives the improved result p^'^\x) 

d(x) 



p^'Hx) = p^'Hx) 1 1 - _ [(fc - -K,,) * (p(i)d)] (x) [ . (6.9) 

The parameter 



K = max 



d{x) 



(fc-lA-2,)*(pWd)](x) 



(6.10) 



gives an estimate for the accuracy of the first order approximation. We will show in Section IVIII that for a small 
quench the entire dependence on the initial state is encoded in p^^\x). Short-range spin-spin correlators are then 
expressed in terms of integrals involving p^^\x) and other known functions, which depend only on the anisotropy 
parameter A. 

A. Example: quench from the Neel state (Ao = +00) to finite, large A 

In this section we consider a quench from the Neel state, for which il-^oci{p) is given in p.2ip . The lowest order 
result for /ci'^'(x) (|6.7p is given by 



(j^-,^ ^ sin (2x) tanh (77) sinh(?7) K.v^ sin (2a;) 

[cosh(77) — cos(2x)] [(cosh(77) — cos(2a;))2 + tanh^(7y) siii^ {2x)]d{x) cosh^(77) 



This is proportional to 1/ and for sufhciently large A indeed small. 



VII. CORRELATION FUNCTIONS 

In the last years there has been tremendous progress in the calculation of equal time correlation functions in the 
spin-1/2 XXZ chair^i^i^i^. These results can be applied also to the generalized Gibbs ensemble of interest here. 
In particular, we may use the explicit expressions for short-distance correlators given in Ref. [l^l, which involve three 
functions (p{p), uj{pi,p2) and uj' {pi, ^2)- Examples are 

K.5)^cth(,). + -^, = (7.1) 
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where uj = Ld{0, 0), ut'^ = dxOj'ix, y)|a:,i,=o- The corresponding expressions for the GGE are identical, but the functions 
^p{^j), uj{fii,fi2), ^'{p-1tI^2) are different. Following the finite temperature case, we define functions 



9] 



Q = 



It follows from Eq. (|4.29p that g'if{x) satisfy the integral equations 



where 



E 



sgn(n)e^ 



{x) - 77 



9f, 



^4 



1 + b 
9i 



l+b" 



{x) 
(x) 



9' 



/+ 



1+h" 



l+b 



(x) 



(x)- 



l+b 

n'+ 



l+b" 



(^) , 



i±ix)=i{x^±irj) , c±ix)=± 



4 cosh [rjn) n^^oc ^ cosh {rjn) 

The three functions entering the expressions for short-range correlators are then given by 

, , p dx/ gZ^^ix) g%{x) 



27rVi+r'(x) l+h-i(x) 



'^(Ail,A'2) = W(o)(/il,/i2) - 



g%^ (x) , (^') 



l+b-i(a;) l+r\x) 



where 



4^1 



1 + b-i 1 + r' 



1 + b-i 



(-4^2) - ?? 



-4^1 



1+b" 



(7.2) 



(7.3) 



(7.4) 



a;(o)(Mi,Ai2) = -4fc(i^i - 1^2) + K2n{ilJ.i - ^^^2) 

W(o)(M1,A<2) = -477^(«Ml -«At2) +??-ft^2(Mi-M2)(*^)- 



(-*A*2) , 
(7.5) 

(7.6) 



A. Perturbative results for "small" quenches 



For small quenches as defined in section I VII we may use the iterative solution of the integral equations for b to 
obtain approximate expressions for the functions (p, uj and a;'. Using the parametrization (j6.2p and the lowest order 
solution (|6.7p . we find the following expressions for the first order approximations to the functions (|7.5p 



W(1)(M1,M2) = W(0)(Ail,/^2) +2 
'^(1)(M1,M2) = UJ[:^){^ll,^l2) +^l 



dx 



— d{x + ifii)d{x + ifi2)p^ '{x) , 

TT 

dx f ~^ dx 

— d{x + 1^12) p^^\x)c{x + i 111) ~ rj / — c(a; + 1/^2)^^^^ (2^)^(2; + i/^i) ,(7.7) 



where we have defined c{x) = c+{x) + c_(x). We can determine the functions entering the expressions for short- 
distance correlation functions such as (j7.ip by numerical integration of the appropriate partial derivatives of (|7.7p at 
(0, 0), once the function is known. 

Results for a quench from the Neel state, where p^^\x) is given by (|6.1ip . are presented in Table ID We stress that 
these are not obtained through a numerical solution of the nonlinear integral equations; we have merely evaluated 
the analytic expressions (j6.1ip(j7.7|) . To the best of our knowledge, these are the first closed- form expressions for 
correlators in the stationary state after a nontrivial interacting quench. 
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A— 2 

(k~0.1) 


i 




(o"io"r+<?> 


A— 3 
(k~0.03) 


I 








1 


-0.661371 


-0.338629 




1 


-0.815293 


-0.27706 




2 


0.376573 


0.056895 




2 


0.643582 


0.039557 




3 


-0.279034 


-0.009872 




3 


-0.57034 


-0.005788 


A=4 
(k~0.01) 


/; 






A=5 
(k~0.008) 


^; 








1 


-0.887611 


-0.224779 




1 


-0.925316 


-0.186711 




2 


0.779648 


0.025859 




2 


0.852509 


0.0177252 




3 


-0.730205 


-0.003019 




3 


-0.81808 


-0.001699 



TABLE I: The short-range correlators for various values of A after a quench from the Neel state (Aq = -f oo) at lowest order 
of the perturbative expansion; k (|6.10p is an estimate of the relative error for p^^^ix). 



VIII. INITIAL STATES WITH NONZERO LONGITUDINAL MAGNETIZATION 



If the initial state has a nonzero magnetization 



(8.1) 



the integral equations (j5.5[) must be modified as follows 

logb(x) - logb(a;) + /i = + k) * log(l + h)](x) - [(fc_ + k) * log(l + l)\{x) 



dx 



k * 
k * 



91 



l+V 
.9; 



[x)- 



1+b 



k+ * 



9u 



1 



{x) 
(x) 



'-^d{x)(^^ + - 

TT ^ '\l+h-\x) 1 



2 dx / g^^jx) 

^ TT Vl+b-i(a:) l+r\x) 



9o (x) 



■2) 



where h is the Lagrange multiplier of the total spin in the longitudinal direction. 

For a small quench, for which the magnetization is close to zero, the right hand side of the last two equations is 
still small, so we can still use the approach developed in Section IVTl The main difference is that at lowest order the 
function (^(x) is now a nonzero constant. The modification is in fact trivial, and we easily obtain 



sinh(C/2) 



I^.^dix)pi'\x) 



cosh(C/2) ' 



(8.3) 



where Pq\x) is the function (|6.7p . As a matter of fact, at the lowest order a;(i)(/xi, /i2) and w|-^j(/xi, /X2) have the same 
form in terms of the generating function £7.1,^ as for = 0. The only difference is in ip{fJ.), which is now different 
from zero 



(/3(i)(//) =tanh(C/2)[d*p^^^](-i/x). 



(8.4) 



IX. SUMMARY AND CONCLUSIONS 



We have considered the problem of the late time behaviour of short-distance spin-spin correlators in the spin-1/2 
Heisenberg XXZ chain after certain quantum quenches. We focussed on the antiferromagnetically ordered regime of 
the XXZ chain and assumed that at infinite times a generalized Gibbs ensemble is reached. Following ReL [4^ we 
constructed quantum transfer matrix description of this ensemble, and derived a set of nonlinear integral equations that 
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describe the largest eigenvalue of the quantum transfer matrix. We emphasize that there is an important difference 
between our description of the stationary state after a quench and thermal ensembles for generalized integrable 
Hamiltonians of the kind considered in Ref. [i^ : in our case the "input data" are not given by the Lagrange multipliers 
(inverse temperature, chemical potential, etc.), but instead fixed implicitly through the expectation values of local 
conservation laws. 

Our main results are as follows. 

(i) We have presented a method for calculating the expectation values of local conservation laws in simple initial 
states (of product or matrix-product form). We obtained explicit analytic expressions for two cases of interest: 
the state with all spins aligned in the transverse direction (Section IIII A[) and the Neel state fSection llll B[) . 

(ii) We have shown that as long as the expectation values of local conservation laws arc known explicitly, it is 
possible to avoid having to determine the Lagrange multipliers defining the density matrix of the generalized 
Gibbs ensemble. This substantially simplifies the calculation of short-distance correlation functions in the 
stationary state (Section |V|. 

(iii) We showed that for "small" quenches, as defined in the main text, analytic expressions for correlation functions 
can be obtained. The properties of the initial state enter through a single function (Section IVII) . which we 
determined analytically for quenches from the Neel state (Section I VI Ap . Using the results of Ref. [47| for finite- 
temperature correlation functions in the XXZ-chain, wc obtained expressions for short distance spin correlation 
functions. In particular, we obtained an analytic expression for short-range correlators after a quantum quench 
from the Neel state to an XXZ Hamiltonian with sufficiently large anisotropy parameter A fSection lVII Ap . The 
next orders of the perturbation expansion in 1/A^ are easily accessible as well. 

Much work remains. For other initial states such as jx, f) considered in section lIlIAi we need to solve the nonlinear 
integral equations numerically or implement an iterative analysis similar to that for the Neel state. Work on this is 
in progress. For quenches to the critical XXZ chain |A| < 1 the role of the conservation laws discovered recently^ 
needs to be clarified. 

While this work was being written up, a preprint by B. Pozsgay appeared on the arXi"\*^, where the same problem is 
studied by means of a quantum transfer matrix formulation, that allows for the calculation of short-range correlators 
by using the corresponding finite-temperature results. The main differences to our work are the treatments of the 
initial state and of the Lagrange multipiers. 
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